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ABSTRACT 



This thesis is a theoretical and experimental study of a 
scheme for convCj.ting energy from a streaming fluid through an 
elastic plate to an electric circuit. The model is a quasi-one- 
dimensional motion of the elastic plate under the influence of 
the supersonic fluid and a transverse magnetic field. 

Anmalysis of the model leads to the equations of motion of 
the plate in the unrestrained case. By assuming traveling wave 
solutions, the dispersion relation, resulting from the coupled 
equations of motion, is solved. The boundary conditions are dis- 
cussed and the resulting eigenvalue problem is stated. By the 
method of coupling modes, the effect of the coupling on the wave 
numbers is determined. 

In the experiment, an elastic media is tested for response 
with electromagnetic and mechanical excitation. The supersonic 
flow is attained with a set of convergent-divergent nozzles 
feeding a test section. The combined plate and flow experiment 
is run and results discussed. 

The conclusions are; (1) the theory predicts an instability 
in the natural modes; (2) the plate can be excited in the second 
modes and the frequency response increases with tension. (3) The 
effects of the instability in the combined experiment are masked 
by forced excitation probably due to the turbulent boundary layer. 

Thesis Supeirvisor; James R. Melcher 

Title; Associate Professor of Electrical Engineering 
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CHAPTER I 



Introduction 

Energy conversion has historically been accomplished by rotat- 
ing machinery, making use of small gap lengths and conduction in good 
conductors. Recent investigations have lead to novel static and 
dynamic converters. The feasibility of generating d.c. povjer from 
the interaction of a moving conducting fluid and an electromagnetic 
system has already been demonstrated. Several schemes for A.C. 
power generation have been presented, however, these are restricted 
by the requirement of large magnetic fields and a good conducting 
fluid. An alternative method for converting energy is shown in 
Figure 1.1. It is a device for producing A.C. power by abstracting 
energy from the fluid flow. 

The device consists of an elastic plate in a channel with a 
fluid flow on both sides. The plate is constructed of highly per- 
meable material and is free to move in the transverse direction. 

The channel^ also of highly permeable material, contains the fluid 
and copper conductors. One set of conductors is excited by a con- 
stant current (D.C.) source, causing a transverse magnetic field. 

The interaction between the fluid and the plate causes transverse 
motion of the plate, similar to the ’’flutter'* motion in aero-elasti- 
city. This motion provides a variable coupling of the first cir- 
cuit with the second magnetic circuit, producing A.C. power, as 
Indicated in the lumped parameter model shown in Figure 1.2. 

Analytically the problem of the thin airfoil performing small 
lateral oscillations in a uniform stream of incompressible fluid 
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was for many years the heart of all flutter prediction. The com- 
plete solution was first published in the United States by Theo- 
( 2 ) 

dorsen . This approach has been further developed by many workers 
and is presented most comprehensively in the text, Aeroelasticity . 
by Bisplinghoff , etc.^^^ 

The definition of panel flutter is the self excited oscilla- 
tion of the external skin of a flight vehicle when exposed to an 
airflow^on one side The analysis of panel flutter has been 
pursued in great detail recently, but there is still some disagree- 
ment between theory and experiment. It appears that the theoreti- 
cal considerations to date do not approach the situation described 
here. The tensions on the panel are significantly less than those 
to be applied in this device. The mechanism for flutter is not 
explained physically in the context of a feedback situation. How- 
ever, the effects of aspect ratio, elastic foundations, damping 
and boundary support have been covered in detail. In the appendix 
a comparision of the prediction of flutter theory and the pre- 
dictions of this theory is shown. The aeroelastician in general 
has confined his analysis to that of "infinite” flow, whereas in 
this scheme the confining effect of the rigid walls will be impor- 
tant. In dense fluids, there has also been investigation of flow 
induced vibrations of fla 
follow? 
recent! 

elastic. The latter is concerned with a statistical repre- 

sentation of the turbulent boundary layer as a forcing function. 

Here a-~ain, the models have Infinite geometry and, therefore, are 
not directly applicable. Lastly there are several "handbook" solu- 
tions available for determining natural modes of vibration, for 
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example, see Den Hertog 

/ON 

Melcher'^ has studied field coupled surface waves for both 
electrohydrodynamic and magnetohydrodynamic systems. Ketterer, 
studying under him, investigated the interactions between a fluid 
stream and a fixed structure capable of supporting traveling 
waves By theory and experiment, it was shown that energy could 
be extracted from the flow in the form of oscillations. 

In this study there is contact between the fluid and the 
elastic structure which was not present in the experiment conducted 
by Ketterer. The effect of this contact, a boundary layer which 
could be turbulent, may be considerable in comparison to the 
effect of an instability in the natural modes. 

There are several methods of exciting the plate available. 
First, it could be located transversely in the flow path, simi- 
lar to holding a piece of grass between the thumbs and blowing. 
However, the length of interaction would be limited and therefore, 
the amount of energy taken from the flow would be limited. 

Secondly, the plate could be fixed at one end only, and the exci- 
tation would be dependent on the trailing vorticics. This would 

(4) 

be similar to the study reported by Eagleson. Another method 
that has been suggested is to place the plate in a tube closed at 
the downstream end and put the combination in the supersonic flow. 
This would provide greater excitation since pressure variations 
considerably above stagnation pressure have been observed in such 
a configuration. The overall efficiency of this scheme is ques- 
tionable, however. 

Tha method chosen was considered both theoretically and experi- 
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mentally the siraolest to study, with the added feature tiiat the 
sis?nif leant Interaction length might produce a reasonably effi- 
cient ’device . In Chaoter II, the device is studied by assuming 
a quasi-one dimensional motion of the olate. The resulting dis- 
persion relationship is solved in chapter III and the implica- 
tions discussed. In chapter IV the longitudinal boundary condi- 
tions are treated. The experiment and results are described in 
Chapter </. 
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CHAPTER tr 



AMLYSIS OF DEVICE 



The device will be modeled by assuming a quasi-one dimen- 

I 

sional motion of the nlate. The model, as seen in cross-section 
from a top view, is shown in Figure 2.1. The flow will be assimied 
Inviscid and ^incompressible in this analys:iS since the objective 
is a device *that can be tested in a wind-t mnel. The proposed 
orientation will also allow neglecting the effects of gravity. 

The fluid flow, the olate and' tlie magnetic effects will be treated 
seoarately, Lnd thea the results will be combined to obtain the 
equations of motion 

2 1 The Fluid Squat iOrtS' !^ " ' 

' ‘ ‘ \ ■ -n-_ ■ I 

In the bulk of the^, fluid, if body forces are neglected, the 



momentum equation is /V 



/ 



> /. 



Du 

Dt 



1 4 . 






- 0 



( 2 . 1 ) 



^ \ f i'- . L > ' : i C 



f- r 



for an inviscid fluid and the equation for continuity of mass is 



^ 4- p v-u. = 0 
Dt 1 



( 2 . 2 ) 






”3 ft 



for comoressible flow we"*wi 11 assume small oerturbations and 
that the waves are long compared to the significant dimensions, 
that is; 
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analytical model 
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v^* “ v' (x, t) - V 

X X ’ o 



V* = V (x,t'> 

y y 



p(x,t) = p' (x,t) “ 



p(x,t) = p'(x,t) - pQ < < 

“T 



where the orimei quantities are the actual values and the unpri-^ 
med quantities the oerturbations.-j-lnserting these into equations 
(2 1) and (2.2) and retaining only fi)rst order terms: 1 



<)v 









V 






f 6 ^x 

t:i 






1 <. 

I ' / i 



Pj 



H ^O ^ ^ P " ^ 



X . Z \ 4. f h. 

'^X ^ V 



(2.3) 



(2.4) 



Initially, the thermodynamic process will be considered adia- 
batic, which allows relating the i;har>S®s oressure and density 
^ ilAf. E iL t , 

by a constant. .. _ « . 



"^p = C~'>o 

”:r ‘^R 

2 . 2 Magnetic Eouations 



A section of the magnetic circuit is shown in Figure 2.2. ® 
We assume a perfectly magnetic material so that the coil will i 
induce a magnetic field in the circuit, with the H field existing 
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MAGMETIC CIRCUIT 



FICiURE 2.Z 
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in the gap only. 



From Ampere's circuital law we get: 





H-d® 



Nl 



or H = 



NI 

d+^ 



Nt 



1+^/ d 



Define: ~ steady state H field. Therefore, the 

magnetic fields with reference to figure (2.1) will be: 



o 

H « H 
y o 



iH/i 



where the superscripts refer to the resnective sides of the plate. 
There are several methods available for computing the force on 
the plate. The most direct method is to use the Maxwell Stress 
Tensor described in detail in reference. '' ' The force density 
in the y direction acting on the plate is the sum of the stresses 
acting in that direction on each of the faces of the plate. The 
stress tensor: 

"Ij * - 2 «lj«kk 



In our case, where is the only significant component, this 
becomes 
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and the force per unit area is 



f « ^ 

y 2 o ^ ^2 



(l+?/d)" (1-^/d) 



Keeping only the linear tenns, the force density is sicrnly; 



■ 4 t 1 , 

f -2n H 5 

y ' o O d — . . 



Therefore, the force is proportional to the displacement and will 
have a destabilizing effect on the plate. 



i 2 3 The Plate Equation 



The Plate has a density p and a modulus of elasticity E. 
Using simple beafli theory ^ the bending moment at any point is pro- 
portional to the rate of change of the slope. 
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the forces and moments acting on it. 
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m = bending moroent. 



Now, surcming the Kioments about the lower left hand corner: 



-dm - vdx + pd4 0 or 



dm , d' 
ox dx 






w , 3 
dx 



dx 



and the forces acting on the element, 



- , dv 

f + = pm2A — ^ or 

ht 

, + p i!i . £i d 

dx^ ” dx^ 



(2.5) 



The three previous developments can be combined to give the equa- 
tions of motion. 



2 . 4 The Equations of Motion 

Equation (2.4) can be rearranged and integrated directly. 
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(2.6) 



The transverse fliaid velocity must be zero at the channel wall 
in both the uoper and lower regions: 
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If one assumes no cavitation the transverse fluid velocity must 
equal the plate velocity at the nlate. 
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the above two equation* maybe added: 
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Equation (2.7), the plate equation, maybe rearranged and 

111 ^ 

inserted in (2.3) and (2.8) to eliminate (p -p ), the pressure 
difference, from our equations*. 
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where we have defined V * 
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It will be convenient to have the equations in nondimensional 
form before further manipulations: 



Variables: , x®x/d, k-kd,o)='a)“ , t t 

o 



P ^ 

^ o 

Parameters: ix - , density ratio, 



H - — , Fluid ^lach number 

a V ’ 

o 

1/2 

Material ”Kach*' number 

m V 

o 

7 

n H"d 

H ^ , Magnetic pressure 

pmv^A 

F T 

D = , Rigidity 

w2prav^d A 



We have normalized the narameters to a dynamic pressure 
1 2 

( 2 P then multiplied by the density ratio. This will 

simplify the solutions as well as provide an indication of the 

strength of the various forces or pressures compared to the fluid 

dynamic pressure. « 

dp^c^ 

V7e now multiply equation (2.9) by — 
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and equation (2.10) 
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( 2 . 12 ) 



The above form the equations of raotion of the system in nondlmen- 
sLonal form 

To solve these equations as they stand for specified forcing 
function would be difficult Indeed. This approach has been fol- 
lowed on simplified equations, where the forcing function was a 
statistical representation of the boundary layer excitation. The 
approach we will follow is to look for instabilities in the free 
inodes, when the source of oscillations is not a lossy one as it 
is in a turbulent boundary layer. 
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CHAPTER III 



THE DISPERSION EQUATION 



The dispersion equation, or the relationship between m and k, 
will give an indication of the stability of the plate in the free 
modes. We will assume a delta function in time and space as the 
excitation, so that the equation of motion will have no driving 
function. In order to determine whether the structure can suo- 
port traveling waves, we also assume that the solutions will be 
of the form.^ _ ^ and.^ „ ^ gj(^*>t-kx). Insert- 

ing these into the equations of motion we get: 



(^>-k)(pK^ - ^ (K^k)v = 0 

a m ' 2 a ^ 



(3.1) 



-k (H-Dk^- M^k + ^ (co-k)^ = 0 

in ^ 



(3.2) 



Setting the determinant of the coefficients equal to zero 
will insure that the equations of motion are compatible: 
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This is the desired dispersion relationship, and its solu- 
tion, 03 ~ f (k) , and the inverse, k * f (eft) , will provide insight 
into the stability of the infinite plate. Before proceeding with 
the solution, there are several limiting cases which should be 
investigated. These in effect will separate the two systems. 

1) If we allow the tension, P, on the plate to become 

infinite, the plate becomes, a rigid boundary. 

Thus, if M — , we get: 

’ m ° 

+ M^k^ (-k^) =0 

or, in dimensional form: 



2 2 2 
(o3-V ic) = M k 
o a 



(3.4) 



which is the expected dispersion equation for convective acoustic 
waves 

2) Similarly if we allow the density of the fluid to approach 
zero, we find the equation of the plate alone. 
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-T 

It is comforting to note that the magnetic force will have a 
destabilizing effect, while the tension and rigidity tend to sta- 
bilize the solution. 



3) Allovjing the speed of sound to become infinite will res- 
trict the solution to the subsonic region. In dimensional form: 




£1 
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+ P 
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2pmA — ^ + 2p V ^ 
^^2 1-0 o d 
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This is the equation for the plate with a dynamic pressure term 
which can be derived from Bernoulli's equation in a linearized 
form (assuming ^ — ► 0) 
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then : - 

- 2p i ■ 

*^o o d — 

The above three checks indicate that the model and the ana- 
lysis are consistant. It also appears that they are to some 
extent physical. - - 

» 

The third check gives a static instability which has been 

observed at subsonic speeds, , 

- - ■S ■ ! . , 

T 

The uncoupled eqxiations (3.4) and (3.5) can be solved simply, 
since they are quadratic, and plotted, as in Figure 3.1. For the 
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rest of the analysis, some experimental considerations will 
govern the approach. Supersonic speeds are difficult to attain, 
and in general, the experimental apparatus is designed for a 
specified flow Therefore, we will take to be a constant 
ec^ual to twice the speed of sound. The rigidity of the plate is 
also a fixed quantity for a given experiment and the effect of 
its variation would be roughly equivalent to varying the tension 
on the plate. Since the plate tension is easily changed experi- 
mentally, it will be allowed to vary in the solution of the dis- 
persion equation. 

We now rearrange the dispersion equation as a polynomial in 

CD •. 



m^-hD^(-2k)+ 0 ? (H-Dk^- + k^ ) 



m 



a a 



+ m (.?k) (-H + Dk^ + M^k^+ M^) 
' m a 



+ (-k^)(M^ - 1)(H-Dk^ - M^k^) » 0 



(3.6) 



and, using the computer, solve for real values of k. 
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The results are shown In Figures 3.2, 3.3, and 3.4 for the 
following cases: • ■ I 

i» 



1) 


Figure 3.2 


^ if - .875 
m 


(strong tension) 


2) 


Figure 3.3 


M = : 545 
m 


(medium tension) 


3) 


Figure 3.4 


M - .314- 
m 


■ (weak tension) 






The solutions plotted^in Figure 3.1 remain as the asympotes 
for these solutions RetT.nrning to equations 3.3 and 3.4, v;e find 
in dimensional form*. 




0 ) ('/+c)k ' I for the fluid 

•- X ■ 

pl/2 ' 

(ju=i (n for the plate 

In the case of strong tension the coupling is weak. Kewaver, as 
the tension is decreased- as in Figure 3.3 to a point where , 

the wave speed in the plate j is approximately equal to (V-c) , 
the "slow* fluid wave, the coupling becomes pronounced. Decreas- 
ing the tension further, decreases the coupling until the solu- 
tions are uncoupled ‘ 

The two systems, the fluid and* the plate, are in themselves 
stable at at least if the majne tic" field is small. The above 
solutions indicate that when coupled together an unstable situa- 
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tlon is produced. That is: od has a negative imaginary part 

, s at+j(cjot-kx) 

and the assumed solution: gj has the form: „ ne 

such that the displacement grows in time. This is the result 
of a feedback mechanism within the coupled system. Let's assume 
a deflection of the plate in the second mode, with a positive 
amplitude over the leading half of the plate. This provides a 
compression wave in the fluid. As this disturbance is swept 
downstream, it arrives over the trailing half of the olate in 
time to reinforce the motion of plate. The^ the motion is trans- 
mitted 




back upstream through the plate, to arrive at the right time to 
start the procedure over again. This is the case when the velo- 
cities / , and V-c are equal. Since this is a distributed pheno- 
m.+ 

roenon over the half length of the plate, the mechanism will work 
for velocities that are approximately equal as well as exactly 
equal. 

A complete discussion of the dispersion equation as it per- 
tains to the infinite problem would not add significantly to the 
study of the device proposed. However, the steps will be outlined. 
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Having solved for a> == f(k -), the coimlex values of cu obtained 

real 

are substituted and the equation solved for + jk^ * 

to determine the type of instability Involved. This can be done 

using the Bers-Briggs criteria , with the help of computer 

technique giving a graphical display of both co versus oo. and 

( 12 ) ^ ^ 

k versus k. as described by Mills. 

The disoersion eouatlons shown in Figures 3.2 - 3.4 have 
characteristics which indicate that the basic theoretical model 
is physically meaningful. The coupling between the flow and the 
elastic media occurs at long wave lengths, with the short wave- 
lengths stable and uncoupled. These are the properties of the 

stream structure previously investigated in the context of an 

(81 

electrohydrodynamic Interaction. There it was shown that the 

effect of the longitudinal boundaries was to make the instability 
absolute. That is, the theory supports the conjecture that grow- 
ing oscillations can be produced. 
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CHAPTER IV ] 

I 

I 

I 

BOUNDARY CX)iTOlTIONS f 



I 

The aoolLcfttioa of the boundary conditions will be aooroachei 
by two methods, j an exact method resulting in the eigenvalue pro- 
blem and an aoproxltnate method called the coupling of modes 
approach. A si^.ple linearized examote will be covered first fol-' 
lowed by a discussion of how to ao'ly the boundary conditions, 
j This leads to the statement of the eigenvalue problem. Next, the 
couoling of mod^s is ■dtg^cussed^. This will give an indication of 
a more meaningful way of apoTylTig*the iSoiTrfdary 'coifditiotis', and 
an estimate of the coupling effect. The boundary value problem 
i is not solved completely. 

T ♦ t- t" “5^ ^ • - F 

I 

O O* f - A 

I 4 1 A Simple Example 

Consider a membrane fixed at the ends. This is similar to 
the elastic plate if the rigidity and the fluid are neglected. i 




Equations of motion: 





2, Assume Solutions of the form: 



^ “ (A sin k +3 cos k ) 

X X 
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A MEMBRAWE FnCEO AT THE ENDS 



FIGURE 4.1 
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4. Soundary Conditions: 

1) ?(0,t) » 0 3 = 0 

2) ^(L,t") = 0 A sin kL = 0 and k * ^ 



5. Solution: 

^(x,t) 



Re 



f 



(An sin 



nirx 

L 




4 2* Method of Characteristics 

tn the example the basic steps to the oroblem solution are 
stated. There is a significant change when the fluid is added, 
especially when it is allowed to have a gross translational 
motion The question arises - at what point may one apply boun- 
dary conditions on the flow such that they are meaningful? V/ere 
the fluid motionless, the boundary conditions must be applied 
at two points in space. This is also true if the velocity of the 
fluid is less than the speed of sound in the fluid. From the dis- 
persion equation, however, the Instabilities occur at supersonic 
speeds only. 

In order to understand this more clearly, the method of 
characteristics is useful. This method uses the general solutions 
of the wave equation for the transient case: 



f = f + (x-ct) + f - (x+ ct) 
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where f+ is a wave -propagating in ths_plus x direction and f- a 
wave in the negative x direction. The solutions are represented 

r . ■ 

by tbtf third dimension of a graph -sphere the first two are time and 

■'» .'/ / // / f, '.'.j//, 

I space. Thus a disturbance’ can’ I5e"followed lij'both time and space « 



Let : 



1 « X + ct and B = x - ct 



Assume a disturbance: p = P 



- a < X < + a 



= 0 elsewhere. 

for the fluid flow shown in Figure 4.,^. 

■ 

In Figure 4.3 the disturbance is shown. 

Lines representing constant and\0 can 'be dravm on the 

\ 

graph Tliose originating at x = + a are>^thG ’most interesting 

since they show the progress of the disturbance. For example, 

^ * 

at t = t the position of the ^disturbance cart be determined from; 
f = f+(B) t = t + f-(^> t = t 

O ' O / . 



The ilope of the lines is of interest since, if the lines of B 

had a oositive slope, the point x -L, would never know the dis- 

\ 

turbance This is the case when the fluid is allowed to have 
motion at a velocity greater than the soeed of sound. For a trans- 
lating media, we must redefine: 

‘'t I'lD wA\^f ■ 



n * X +(v+c)t 



and B = x+(v-c)t 
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as the characteristic solutions to the convecting wave equation. 
Then (v-c) , the slope of the B lines^ becomes positive for v > c, 
and points upstream will not be effected by events downstream. 

This implies, in the scheme prooosed, that the fluid boun- 
dary conditions must be applied at the leading edge of the ulate 
in order to determine their effect on the interaction. The fluid 
will not know what is happening at the trailing edge until it 
gets there. There are tv;o variables on the flow, the pressure 
and velocity, and their values at the ooint x * 0, as a function 
of time, form the constraints on the flow. 

4 3 The Eigenva lue Problem 

The dispersion equations was sixth order in (u . The six 
boundary conditions are the two snecifled for the flow, plus the 
four on the plate From simple beam theory the latter are deter- 
mined from the tyoe of support- free , pinned, or clamped. If we 
assume the ends are clamped, this implies that that the deflec- 
tlon and slooe, , are zero at x * 0 and x - L. Since v(x,t) 
and p(x,t) are defined as difference quant5.ties (above and below 
the plate) , their values are zero at x = 0. 

Returning to the equations of motion this can be stated sim- 
ply. 

(m - k)(|j>la-P.E.)q ‘1 I V « 0 (4.1) 



k(P.E.)n + I (m-k) V - 0 



(4.2) 
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(P.E.)ri + p - 0 . (4.3) 

where P. E. * H-Dk'^-M^k^ + 

m 

As in the example, assume a solution of the form: 

T)(x, t) = He ^ A- e j 

and rewrite equation (4.1) and (4.2). 

P = - ^ n (4.4) 

V - ^ ^ . f „ (4,5) 

k(P.E.) ^ 

apo lying the boundary conditions; 



(1) n(0,c) - 0 ^ - P (4.6a) 

(2) ^ (0,t) - 0 f ■■ 0 (4.6b) 

(3) T^'(L,t) “0 ^ ^ e^^n^ « 0 (4.6c) 

(4) ^ q(^,t) =* 0 ^2jk A(i e^^n^ 0 (4.6d) 

<JC ^ 



^ ■ (P.E!(k»,) ^ 



(5) p(0,t) « 0 



(4. 6e) 
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(6) v(0,t) =0 £2f, (k ) A = 0 (4.6f) 

in n 

Rearranging in iratrix form; 

A(o),k)* 



1 


1 


1 


1 


1 


1 




^2 






^5 


^6 








e^V 








k2e^^2^ 
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k/e-* 4 
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k3e^^5^ 


keeJ^^e" 
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1 


■ P.E (Uj) 


■ P.E.(k'2) 


■ P.E.(kj) 


■ P.E. (If^) 


■ P.E.OTj) 


■ P.E. (Eg) 


£,(ki) 


f.j^(k2^ 






fjCk,) 


fj(kg) 



In theory then the equations = 0 and D(oa,k) »= 0 can 

be solved. Assuming a complex oa, the six k's are found from 
D(co,k) * 0. These values are substituted in to A, This procedure 
is continued until A =0; both it's real and imaginary parts, 
since it is complex In practice another method, that of weak 
coupling of modes will be used to attack the problem. 
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4.4 ^^ethod of Coupling of Modes 

One might ask the question: *’lf two systems are individually 

stable, why would a weak coupling between* the two result in totally 
different motions?*' The answer is that the basic solutions for 
the two systems are approximately correct, and that the changes can 
be predicted within certain limitations. First, the systems must 
be stable by themselves. In this problem, this restriction will 
not allov including the effect of the magnetic stress on the 
plate since it has destabilizing influence. Secondly, the coup- 
ling between the two systems must be weak. Within these restric- 
tions, mutual coupling of modes will allow simplifications based 
on a physical understanding of the system. 

(13) 

In theory, the approach is quite simple. A complicated 

system is first divided into its component parts. The equations 
of motion for the isolated systems are solved exactly, and the 
solutions are expressed in the ''normal modes' of the element. 

The complex system is then asstuaed to be made up of these ele- 
ments weakly coupled. The motions within the system are described 
by the elemental motions which are slightly perturbed by the 
coupling. 

As an example, in the system under consideration the ele- 
mental parts are the fluid and the plate. Assuming the fluid 
motionless and the plate of negligible rigidity (a membrane) , 
each will propagate disturbances as waves, both a forward wave 
and a backward wave in each system. As seen from the solution 



«c 



Is. ioUfiisBS. 1ft ittOm £jf 

Tffiw <iTf ^ w43 ifu ^ *(^9 

IlI J/Vft#? 0*t7 I^PMI*** « ^t»C«r T«iw 

3«i i ic^4 9uio* iU*» ti^ ^ «p»ciit mfX *^#00^1^ ii»«f If 1 i| 

if«:> r«lft«r«* ftVA 4M^3S|9 9d 

Tt^is f^«am f*)'? ,3 €iI' 9 . fiTiT«> tfk^iftir IrSda^ttvsif »d | 

ril4 •.A4K9 rJ a*#! »aMritl •X4ft9« wH f 

tyii 10^ tiC95^ f •7«fi9i9(«^ 9««9 #i«X ra 7f*> • 

•iiMos* »i:r ,Y£aia»t)»** ^^jD^liau joJi/XXiwwl «(id 3l moIi tatlQ * 

• ^ ,' • 

crvtb < 3 iif 3 iw 2«e*' p« !«»<• %u»i*v4 «4t | j— i»39d |riTf 

l>aod 4i>«l3»lxiX<|<i • iii«r t»ir* ^<9 »i|)f(*«Did f«lrtt<^ «trriC9 ' 

mt«¥« ••fti i« ^ 

I 

^ %««»}• •! «t3 itrroaitl oi 

Mit •71a#' jr®iiWoi!fC»> «9 • »•♦»►• 7ni1 9( 

9*^ L^W Ct9<»l^ *-»» *«9ftl9iai *<19 <Mhi30* lo 

•WMneiia aWt 9gfc— IvvT^e fti f^ t rn gtj tm 9rrt 

I 

-*ife oioji iid «) iHpA«» <i«dj •> aB^an(4 «dJv« mSj 

bM*l9oQAi C7W »<ia 9(fi^ b»fm/<99 flslaM' uSumm 

•Ha ftlPllI* r*I^«i9itt •ifl 

-af« loia MtjauaftitMxa 3^«iii Ajayt »«n tii n* %A 

? ' T • 

>Mn «ff5 99«fi| Ms n<u mTM Oic**! I«3(ia0 

I ^ aftcafi^i^a %) *5Mv >f«3 &;:*• 

9*f<^ k*T*n09<e^ • #994 , f««w< t4 aoKrnrtf^^** ipf# lfi» <fa*9 

0023 ^ foj 94i3 9ft0« tA 4iMa o »»*** a tai i tf i Ins 



39 



of the dispersion relationship, the coupling will occur between 
the forward membrane wave and the ’’backward' fluid wave. There- 
fore, the other two waves can be neglected. The motion of the 
system will be that of these two waves slightly altered by the 
coupling. 

4. 5 Elemental Equations in the Coupled Mode Form 
1,. Fluid Equations 

The fluid equations can be greatly simplified to get an 
approximation to the solution. 

?^v(x.t) ^ _ 1 ^p(x,t) 

f^x p (v+c) ^t (4,8a) 

r Q — 



^p(kvt) ^ _ ^. v . C x^ t ). 

>x Po ht 



(4.3b) 



These coupled equations are most easily solved if a method of 

(13) 

decoupling them can be found Following, Louisell,' the 
decoupled form is derived from linear combinations of the above 
equations. That is; assi.tme a normal mode form: 
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1 2 1 

b ± (x,t) = 2 (P) P(x,t)) 
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where the mode amplitudes are normalized so that the square of 
the amplitudes represents energy carried by the mode. 

Multiplying equation (4 8b) by the fluid impedance ( — ^ 

^o ~ 

and adding and subtracting from equation (4.8a) will give the 
desired decoupled form: 

rQ — — ^ 0 

or 

^ ( b + (x,t) ) ' - Ij (b i (x.f) ) 



Assume that the time dependence will be of the form 
then: b + (x,t) * a + (x) e^^^ + a*+(x) e 

Now if v(x,t) = Re V(x)e'^^*' 

1 __ 

then a ± (x) = (p)2(V(x) + ^ 

and equations (4.8a) and (4.8b) become 
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( ^ “ jk)a? (x) * 0 
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( - jk)a-(x) = 0 



(^.^ •+• jk)a-(x) “ 0 



(4.9b) 
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where k == 

v+c 



2 The P^ate Equations 

Again simplifying assumptions must be made. Neglect- 
ing the rigidity and the external forces on the plate, the nor- 
mal mode form of the equations for this element can be derived 
in an analogous manner. 
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there fore , 



( -jk 2 >a 2 +(x) =0 ( ^ a* + ■= 0 



(4.10a) 



( ^ +jk2> a2-(x) =0 ( ^ -jk2>a* - = 0 
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where 




0 ) 



4 • 6 Coupling Between the Two Elements 

In section B, we derived the normal mode form of the element 
equations In effect the single variables in equations (4.8a) and 
(4 8b) have been replaced by a combination of the two variables. 

To put the previous equations in coupled mode form, it is nec- 
essary to make linear combinations of the elemental forms. This 
is done in such a manner that when the coupling is reduced to 
zero, the equations reduce to those derived for each element 
separately. 
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1 1 / 2 — 

32 = 2 ^ p S 

and 

^11 ^ ^ ^-/4p^k2> 

(, = -i 

32 ^ 

0 A coupling coefficient 

In order to make the problem more tractable we must make use of 

our previous results. The first is that coupling between the two 

systems will occur when the two wave members, k., and k 2 , are 

approximately equal. Since k^ = and k 2 *^ and in Chapter II 

m 

we saw that the coupling occured when *** “ V-c, we can there- 
fore neglect the fast fluid wave (V+c) and the backward plate 

wave (V ) . 
rn~ 

Secondly, we must determine if, in fact, the coupling between 
the two systems is weak. Returning to the coupled dispersion 

relationship and inserting our previo".c assumption that the 
rigidity and the magnetic field are zero, we have: 
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or in dimensional form and rearranged. 






where the term on the right represents the coupling between the 
two systems. Now let the perturbation of the elemental wave 
numbers be small: 
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where 6 , ^ 7 ; — 

1 V-c 



k 2 = + *2 where d 2 <^ 

and keeping second order terms in 6 we find: 
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If we now let 6 . 



0 and recall that v-c *= c, we see that: 



1. The coupling is proportional to ratio of the densities 
of the two media |x, which is a small quantity, and 
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Therefore, = jk 2 and equations 4.11 reduce to: 






B. 
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(4.12a) 
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(4.12b) 



These equations can be solved by assuming solutions of the fom 

X 

e and setting the determinant of the coefficients equal to 
zero . 
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where k- = ^ — and k« ® ^ , the undisturbed elemental wave 
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4. 7 Summary 

The eigenvalue problem has been stated, but was not solved 
because the method of coupling of modes gave better insight into 
this problem. We applied coupling of modes to a simplified pro- 
blem and determined the perturbation of the wave nembers based 
on the effect of the coupling. In order to properly apply the 
boundary conditions in this case, the approach of coupling of 
modes could be followed. But it is not clear that the neglected 
waves can in fact, be neglected. Since at long wave lengths, the 
boundaries will have a definite effect on both modes, this must 
be accounted for in the theory. The results of Chapter III and 
IV do provide a basis for experimentation. If we assume that 
the wave numbers are only slightly perturbed; we can compute 
from the physical length, from equation 4.13 and and 0 )^^ 
from the dispersion relationship. Although these values are not 
exactly self-consistent, they should be approximately so. 
Experimental techniques to achieve these values are described in 
the next section. 
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CHAPTER V 

EXPERIMENTAL PROCEDURE 



In order to experimentally test the predictions of the 
theory developed, it is necessary first to insure that the 
’’ingredients” required are available. The initial attempt at 
assembling a device and making test runs was not successful. 

The step by step, approach will be described in detail. This con- 
sisted o^ several separate experiments. The first was testing 
of the plate structure to be used. The second was investigation 
of a means of attaining supersonic flow. Next, the two were com- 
bined to test the flutter theory. 

f 

5, 1 Testing of the Plate 

There were two experiments performed on the plate along. The 
first was an investigation of various types of material and vari- 
ous end conditions The, test set up is shown schematically in 
Figure 5.1 Due to the magnetic coupling between the two cir- 
cuits, the results of this test were qualitative rather than 
quantitative. The results are shown in table 5.1. The end 
conditions have a noticeable effect on the response. The damp- 
ing or losses ?ff^condition but tightly clamped appeared to nulify 
any response. 

The frequency did increase with tension which indicates that 
the use of membrane analysis in the coupling of modes section is 
not totally in error. The frequency of response was reasonably 
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close to that predicted. (See example 5.1) 
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The second experiment on the plate above was to deterniine 
if in fact the second harmonic frequency could be excited when the 
media had the perdomlnate properties of a membrane. Equipment 
similar to that shown in Figure 5.1, was used for this part of 
the experimentation To avoid the magnetic coupling, one of the 
electromagnets 4as replaced by a speaker for mechanical excita- 
tion with weak coupling. The plate thickness" was reduced to .001 

i- 

inch, to more closely .approximate a membrane . d The equipment is 
shown pictorially in Figure 5. 



In this configuration the plate frequencies were not detec- 
table The membrane frequencies gave excellent response and are 
tabulated in Table 5.3 The variation between predicted and 
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actual frequencies can be explained. First, the actual para- 
meters for the material are difficult to determine exactly. 
Secondly, at the higher values of tension, there is a possibility 
that the grips or end clamps slipped slightly. 

The main conclusion to be drawn here is that the second 
harmonic can definitely be observed. The quality of the response 
decrease with increasing frequency since the higher modes are 
more difficult to exert. Tlie frequency expected when operating 
the wind tunnel could be reasonably predicted from this type of 
experiment 

The responses at 355 and 720 cps are probably due to weak 
coupling with the speaker or other portion of the apparatus, since 
they are independent of tension. 

5', 2 Experiment to Obta in Supersonic Flow . 

The physics of supersonic flow is a very complicated sub- 
ject and is extensively covered in the literature, notably in 
reference 15. Only those principles xd:iich are germane to this 
experiment will be touched on. In general, supersonic flow can 
not be attained directly from an air compressor. It is necessary 
to modify the compressor output by passing it through a set of 
nozzles. For this experiment the nozzles were of the convergent- 
divergent type. The design procedure is to attain Mach 1 flow at 
the throat of the test section and then expand to the velocity 
desired. The calculations for this are shown in example 5.2, 
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based on the set of nozzles shown in Figure 5.^. 
EXAMPLE 5.2 (Reference 15, Table Bill) 



Inlet Conditions: 
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Secondly, the flow in the test section will be greatly 
effected by the build up of the boundary layer. This has the 
effect of a convergent nozzle, which tends to reduce to flow 
velocity toward Mach 1. In the case of Mach 2 flow, the velo- 
city will be reduced to Mach 1 in a length of 31 diameters. In 
this experiment, the diameter or depth of the channel will be 
taken as one centimeter to allow sufficient magnetic coupling. 
Therefore, the length must be less than 31 centimeters; an ini- 
tial length of 20 centimeters was used. (From this length, the 
thickness of the plate for a desired natural frequency can be 
computed as in example 5.1). Based on these calculations, the 
experiment shown schematically in Figure 5.3 a. i was constructed. 

The pressure readings were taken with both the plate in 
position and removed. The readings were taken with a guage and 
compared with a monometer for accuracy. The readings obtained 
are shown in Figure 5.4. 

These results are at a variance with those predicted in 
example 5.2, but they do indicate supersonic flow. The variance 
is due to several factors. First, the ideal area ratios are 
significantly diminished by the boundary layer to an "effective” 
area ratio. Secondly, the capacity of the air ejectors was 
insufficient for running the experiment as desired. Also, the 
pressure measurements were made at the wall of the channel, pro- 
bably in the boundary layer, and not in the center of the flow. 

The latter problem is difficult to circumvent, and in previous 
experiments (reference 15, page 135) good agreement between actual 
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PIQURE 5.3 
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wall pressure and theoretical mid-stream pressures has been 
attained. Using the lower pressures recorded, the speed of the 
flow can be approximated. Based on the readings at 35 psia inlet, 
the flow is Mach 1.47. 

5.3 An Experiment to Test the Flutter Theoiry 

We have determined that the ingredients of the theory are 

attainable. The flow can be niade supersonic, the plate can be 

excited in the second mode, and the frequency of this mode varies 

with the tension. Thus the results of Chapter IV indicating 

that the perturbation on the wave number k can be estimated from 

2 

the coupling seem reasonable. The rigidity term Dk can be 

2 

neglected compared to the membrane term at the frequencies in 
question. This is especially true if the very thin (.001") plate 
is used. The first experiments with the combined ingredients 
used a .015 inch plate. The equipment consisted of the nozzle, 
a plate, a magnetic detection device and an oscilloscope. The 
presentation on the scope was confused, showing standing waves 
and traveling waves and these responses were at several fre- 
quencies. These were probably excited by the boundary layer. 

The tension was varied and the frequencies of the response 
increased with increasing tension. By varying the pressure the 
response would set in at a certain pressure for a given tension - 
quantitatively a higher tension required a higher pressure. At 
a given pressure, providing it was above the threshold, the vary- 
ing of the tension was limited mechanically. That is, the tension 
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could be increased until the plate failed in tension, and if the 
tension was decreased the plate would fail by excessive dis- 
placement. The latter occured at low values of tension. It was 
felt that nothing of substance to support the theory could be 
gained from further experimentation with this set up. ’’Flutter" 
did occur and this did provide a variable coupling between the 
two magnetic circuits and there was an output signal. 

However, this flutter could not be attributed exclusively 
to the mechanism proposed in the theory. It was probably caused 
by excitation from the boundary layer. There was a component in 
the vicinity of the frequency expected as computed below: 

(1) length of membrane; k ® “ •628 in 

(M^-l)kv^ 

(2) from dispersion equation: f * = 1080 cps 

-1 

0 )^ = 495 sec 

c3 

(3) from coupling of modes: n * x “ .25 in 

4d^(v-c)^ 



In order to prevent the end connection* from interferring 
with the flow they were removed from the test section as shown 
in the figure. Therefore, the flow is in transition over a short 
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section at both ends of the plate. This could also provide an 
excitation which was not predicted and would be difficult to 
analyze. 

Therefore, another experiment was designed. Due to a time 
problem, It may not be completed before the due date of this 
thesis. The changes Include using a .001 Inch plate to more 
closely approximate membrane performance, refined end connectors 
to Insure a good response and no flow disturbance by the tension 
mechanism. The results of this phase of the experimentation will 
be reported In Appendix B. 
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CHAPTER VI 

CONCLUSIONS AND RECOMMENDATIONS 



6 . 1 Conclusions 

The analysis of the model does predict the desired pheno- 
menon. The Interaction between the flow and the elastic media 
at long wavelengths will theoretically produce grov^ing oscilla- 
tions. The effects of the longitudinal boundary conditions can 
be approached by the coupling of modes. In the simplified pro- 
blem considered in Chapter IV, the effect of the coupling on the 
wave ntimber was shown. This approach could be expanded to show 
the influence of the boundaries on the v;ave numbers. 

The experimental portion shows that the separate systems 
can be achieved. The plate can be excited in the second mode. 

The natural frequencies of the plate increase with tension and 
the tension required by the theory is physically realizable. A 
supersonic flow over the plate can also be attained in the labora- 
tory. With the equipment used, the maximum velocity was 1.47 
times the speed of sound. The major source of difficulty, 
^countered in combining the systems was the excitation provided 
by the boundary layer. This was not included in the model or in 
the analysis. The excitation gave forced responses, mainly in 
the first mode, although there was a second mode in the response. 

There are several other considerations which have not been 
covered. With respect to the elastic media , there will be eddy 
currents induced \dien it oscillates in the magnetic field. 
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These represent a loss in the system \ohlch could be significant. 

The Curie temperature will impose a thermal limitation, restrict- 
ing efficient operation of the device to relatively low tempera- 
tures. From the stand point of strength, it is fortunate that 
ferrous materials, which are required magnetically, also show 
good fatigue properties. 

The most severe limitation on the device efficiency will be 
the losses incurred by the fluid in tr ansiting bhe nozzle and 
interaction area. It was noted in the experimental section that 
the flow would be supersonic only for a short distance. For 
instance, the Mach 1.47 flow in the experiment would be "choked” 
to Mach 1.0 at the end of the test section (six inches), due to 
boundary layer build up. At the same time the isentropic stag- 
nation pressure would be reduced to less than one half its 
initial values. This represents a considerable loss in efficiency. 
A more refined design would improve on this situation, but there 
will be a definte limiting efficiency attainable. 

The subject of sxipersonic flow in short ducts is somewhat 

(15) 

complex and experimental data has not been consistant. 

In general, the boundary layer thickness increases with lengtit. 

But in some cases, negative friction coefficients have been 
attained, indicating a decreasing boundary layer. The effects 
of the oscillating plate on the boundary layer have also been 
studied but there appears to be no significant change in boundary 
layer tlickness. 
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The device itself raay not be the most efficient way of 
obtaining energy from the flow. However, it is possible that 
to employ it in with a conventional energy con- 

verter would improve overall efficiency. The important conclu- 
sion here is that a mechanism is available for transfer of energy 
between the flow and the media. The mechanism itself is not a 
’’lossy" one - the inefficiencies are incurred In the scheme 
for using the mechanism attempted here. In other words, this 
paper has been limited in scope to a particular device, but 
further research would lead to other methods of using this con- 
cept for energy transfer. 

6 . 2 Recommendations 

The study of the longitudinal boundary conditions in this 
scheme could be expanded using the coupling of modes approach. 

The coupling coefficients considered here pertained only to the 
coupling between two waves, one in each system. The coupling 
between the other two waves , both with the other wave in its o\m 
system and the other wave in the other system were neglected. 

The latter still appears reasonable. But the boundary conditions 
will provide coupling between the two waves in one system. We 
neglected this here because at the wavelengths in question, the 
rigidity term could be neglected and the plate considered as a 
membrane. In the general case this would not be so, and the area 
of longitudinal boundary conditions in coupled systems has many 
possibilities for study. 
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The combined experiment showed the o-vowls tending effect of 
the boundary layer, essentially nasking the desired response 
if present. The boundary layer could be reduced either by 
thermal methods or by expanding the flow in the depth dimension. 
If this were done, the correlation of experiment and theory 
might be attained and the mechanism studied. 

(9 ) 

Ketterer used a fluid ._8.tream and a spring coupled 
electrically to study this mechanism. Fraize^^^^ studied a 
device in which a fin ^ situated in a wind tunnel (sub sonic ) , 
vjas excited mechanically to produce traveling waves. Signifi- 
cant reductions in drag were recorded. Replacing the mechanical 
excitation by magnetic traction along the length, using the cot^p- 
ling mechanism investigated here, provides the possibility for 
p\amping of or propulsion through a fluid surrounding the fin 

This would be the inverse of the problem studied here, but 
would not have the boundary layer problem inherent in super 
sonic flow. 
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Appendix A 

A COMPARISON WITH FLUTTER THEORY 



Flutter is ■’defined as the self~excited oscillations of a 
airfoil in sxapersonic flow. This definition leaves some 

doubt as to the physics involved, but since it is restricted to 
supersonic flow it seems conceivable that the mechanism involved 
is similar to the feedback situation described in Chapter III. 

The problem of flutter is stated as a partial differential 
equation, where structural, inertial, and also dynamic operators 
are combined to determine the equation of motion for the panel. 
Two approaches are followed from here depending on the geometry 
involved. One is to apply the boundary conditions to the plate 
only. Since these are two at each end, a "fourth order" eigen- 
value problem results. 

The second approach is to assume traveling wave solutions, 

gjkx-cot^ and ass^mllng real values for k, compute complex values 

( 19 ) 

of a>. Dugyidji^ ’ has used both approaches and compared the 
results. lichen recommends which approach to use based on the 
geometry of the problem considered. 

In terms of our parameters, the equation of motion is; 
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n ^ 



. p A + 2pA + G II 



2 „ 

'’l-a (2i 1. £>i 

a^ a 



(1 



where the term on the right is the 'kerodynamic pressure'.' For 
small M , this reduces to 

3 



2 f ^ ^ V 

■p' 



and the equation is the same as equation (2.10) except for the 

stiructurai damping term G . In our case, the aspect ratio is 

s 

zero and the first approach is recommended. That is: the bound- 

ary conditions on the plate are applied and a solution of the 
form 
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is assumed. This forms the eigenvalue problem, and the solutions 
for CD^ ^ co^ are shown in reference 18. Many parameters have been 
introduced and this leads to lengthy computations to compare the 
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the results of the two theories. For the tension, rigidity, and 
plate thickness employed in the device, the solutions predict = 100 
and 0 )^ = 1000. The approximation is necessary because the 
plotted data must be extrapolated to include the values of ten- 
sion specified. These results show correlation within an order 
of magnitude of those predicted in Chapter V. ( o>^ «= , 

) . Most of flutter investigation has been restricted 
to determining the onset of flutter and methods for avoiding ^^’s 
occurance. Again using our parameters in figure 10 of reference 
18, flutter is predicted to occur when the tension corresponds to 
making ® v-c. This is surprisingly consistent with the pre- 

vious results. The predicted flutter frequency under these condi- 
tions is on the order of 20 cd^ where equals the natural plate 
frequency neglecting tension. The previous results predict 

0)^ =“ lOcu . 
f o 

There is considerable agreement between the two results, 
in spite of the difference in boundary conditions. 
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Appendix B 



The previously described experiment was constructed. How- 
ever, the improved end conditions to lessen the flow disturbance 
were mechanically unweildly. They did not provide the firm con- 
nection desired. Therefore, no further results of significance 
can be reported. 
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